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Abstract 

We analyse the role of the Extended Clifford group in classifying the spectra 
of phase point operators within the framework laid out by Gibbons et al for 
setting up Wigner distributions on discrete phase spaces based on finite fields. 
To do so we regard the set of all the discrete phase spaces as a symplectic 
vector space over the finite field. Auxiliary results include a derivation of the 
conjugacy classes of ESL(2,FAr). 
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1. Introduction 

The Wigner distribution, introduced by Wigner in 1932 [T] for the purpose of 
associating a phase space distribution to quantum systems with as the configu- 
ration space, has, over the years, been extended and generahsed in many directions 
[2]- [5], differing in the aspects of the orginal definition that one chooses to retain. 
Among these efforts, those aimed at setting up Wigner distributions for finite state 
systems have, of late, been of great interest because of their potential relevance 
to quantum information theory 0. For the case of quantum systems with Hilbert 
space dimension N equal to power of a prime. Gibbons et al [5] developed a for- 
mulation for Wigner distributions which elegantly ties up Wigner distributions or 
equivalently the phase point operators with mutually unbiased bases [5, I5| avail- 
able in such dimensions. This approach yields not one but many possible definitions 
of phase point operators and hence of Wigner distributions and Gibbons et al pro- 
pose a scheme for classifying them on the basis of the action of the Clifford group. 
In the present work we examine this classification, taking the eigenvalue spectra 
of the phase point operators as the sole criterion. We highlight the role which the 
extended Clifford group plays in this context. In other words, in line with the 
observations in 9 in the context of SIC POVMS fTUI, we argue that, in so far as 
the spectra of the phase point operators is concerned, what is relevant is not the 
Clifford Group, but the extended Clifford group. This has the consequence that 
the number of 'similarity' classes in the sense of spectra is smaller than that given 
by the Clifford group action alone — for = 7 it is 210 as opposed to 360. 

A brief summary of this work is as follows: In Section II we recapitulate the 
ideas developed by Gibbons et al leading to the definition of generalised Wigner 
distributions. We restrict ourselves to N an odd prime, and set up the notation 
and terminology used later. In Section III, we develop a convenient scheme for 
labelling the phase point operators by elements of Vn+iCZin), an A'^-l- 1 dimensional 
vector space over the finite field Zat, and transcribe the actions of the Heisenberg- 
Weyl (H-W) group and the Clifford group in terms of linear transformations on 
Vn+iCZn). By a suitable choice of basis in Vn+iC^n), the scheme developed here 
is shown to yield a convenient way of labelling the 'affine planes', sets of phase 
point operators related to each other by the H-W group action, as elements of 
Vn-i{Zn). This machinery is put to use in Section IV to study the action of the 
Clifford group on the set of phase point operators and affine planes with the purpose 
of classifying them into orbits under the Clifford group. In Section V, we examine 
the spectra of the phase point operators belonging to distinct orbits and find that 
several orbits are degenerate in the spectral sense, suggesting the presence of an 
extra symmetry. In Section VI we show that the extra symmetry is related to the 
operation of complex conjugation, and that the number of distinct spectra of phase 
point operators can be quantitatively understood in terms of the Extended Clifford 
Group - the Clifford group augmented by the operation of complex conjugation. 
This section contains our principal results. We conclude in Section VI with a 
summary and an outlook. The appendix contains a derivation of the conjugacy 
classes of the group ESL(2,FAr), where Fjv is a finite field. 



2. Generalised Wigner distributions 

We begin by recalling the salient features of the framework laid out by Gibbons 
et al for setting up discrete phase spaces and defining Wigner distributions thereon 
though not necessarily in the same logical sequence. Their general framework ap- 
plies to quantum systems of prime power dimensions, but we confine ourselves to 
the odd prime case. 



1. Hubert Space: Consider the Hilbert space Hn associated with an 
N-level quantum system and introduce in it an orthonormal basis 
{\q): q GZn}, referred to as the coordinate or the standard basis, 
and a momentum basis {\p): p € Zjv} related to it by a finite 
Fourier transformation: 

2. Phase Space: With the iV-level quantum system wc associate a 
phase space consisting of points {q,p) with q,p € Z^, for conve- 
nience arranged in the usual Cartesian fashion. The fact that Zjv 
is a field has the consequence that the phase space can be decom- 
posed into N parallel lines, striations, in iV -|- 1 ways. For each 
striation, we can select the line through the origin as its represen- 
tative and call it a 'ray'. We thus have A'' + 1 rays where each 
ray consists of points of the form {sq, sp) with q,p fixed, not both 
zero and s taking all values in Zjv- For obvious reasons, the two 
special rays (s, 0) and (0, s) are denoted as horizontal and verti- 
cal rays respectively and the corresponding striations as horizontal 
and vertical striations. We arrange the striations with these two 
as the extremes: 

vertical, , horizontal. 

3. Transformations on Phase Space: Two kinds of linear operations 
on the phase space which map the phase space to itself naturally 
arise: 

• Translations, T{qo,po), N"^ in number, have the effect of mov- 
ing a phase point by an amount {qo,Po)- These operations 
constitute an abelian group and map a line in one striation 
to another in the same striation. Thus, for instance, starting 
from a ray, except the horizontal, all other members of the 
corresponding striation can be reached through r(l,0). For 
the horizontal ray, one can achieve this through r(0, 1). 

• 'Symplectic Rotations' 

(^0^(7 0(0' 

with a, 0,'y,6 G Zjv and a5—(3^ = 1 (mod N). Such 2x2 ma- 
trices S constitute the group SL(2, Zjv), containing N{N'^ — 1) 
elements. Its generators are 

^1 ^ ( ' ^' = {-1 0) ■ 

These transformations map rays onto rays, hence the figu- 
rative appellation 'rotations'. 'Symplectic' means that when 
applied to {qi,Pi),{q2,P2) they leave the symplectic product 
{q2Pi - unchanged. 

The semi-direct product SL(2,ZAr) k T{qo,po) will be of interest 

here. 



4- Heisenberg-Weyl group: Returning to the Hilbert space we define 
the operators X and Z on Tijv through their actions on the coor- 
dinate basis \q) (momentum basis \p)) as : 

X\q) = \q + 1) {X\p)^u;-r'\p)), 

Z\q)=u;'^\q) { Z\p) = \p + 1) ). 

Consider the set of N'^ operators 

D{q,p) = t'^pX'^ZP , 

where r = —e^'^l^. These operators, cailed displacement opera- 
tors, obey 

i?(q,p)t = B{-q,-p), 

D{q,p)D{q\p') = T'^^'P~^P'^D[q + q',p + p'), 

Tr[D{q,p)] = iV(5,,o<5p,o, 

Tr[D\q,p)D{q' ,p')] - iV5,,,,<5p,p. . 

This set of unitary operators forms an orthonormal basis in the 
space of all complex operators on T-Ln and furnishes a projective 
representation of the group of translations on the phase space. In 
the standard basis, they can be represented by matrices as 

[D{q,p)U^T'iPuj''P5i,k+, . 

In particular, 

[D{Q, = Z = <5,,fe; [D{\, 0)]^ = ^ = 5,^k+i ■ 

Further {e*'^Z)(g,p)} constitute the elements of the Weyl-Heisenberg 
group (or the PauH Group). 

5. The Clifford group: The Clifford group consists of all unitary ma- 
trices U such that 

UD{q,p)U^ = e'^-^-D{f{q,p),g{q,p)), . 

for suitable phases e^' p and functions /, g It is a semi-direct prod- 
uct of the Weyl-Heisenberg group (playing the role of translations) 
with the group of N dimensional unitary matrices U{F) labelled 
by elements 

F = ^ " ^ ^ ; detF = 1 mod N, 

of SL(2,ZAr). The unitaries U{F), defined up to a phase, act on 
D{q,p) as follows: 

U{F)D{q,p)U\F) = D{aq + 0p,jq + 6p), 

and, up to a phase, are explicitly given by [11] 

^(p. ^ f ^ r^-'("'^'-^^'^+'^')\j){k\ /3 ^ 0, 

1Ej^"^^1«j)(j| 13 = 0. 

We shall denote the unitaries corresponding to gi and g2 by V and 
W. They generate all the C/'s. 

6. Mutually Unbiased Bases (MUB): From the properties of the dis- 
placement operators listed above it is evident that two such oper- 
ators D(q,p), D{q' ,p') commute with each other if and only if the 
'symplectic product' {q'p— qp') of {q,p) and {q' ,p') vanishes. This 
in turn implies that the operators D(q,p) with {q,p) lying on a 



ray commute with each other and hence their simuhaneous eigen- 
vectors furnish an orthonormal basis. Moreover, the orthonormal 
bases corresponding to different rays can be shown to be mutually 
unbiased with respect to each other Since there are N + 1 
rays we are naturally led to the full set of + 1 mutually unbiased 
bases. This association between rays and MUBs is a key ingredient 
in the construction of generalised Wigner distributions as we shall 
see later. As to explicit expressions for the MUBs there are many 
choices depending on the choice of phases for the vectors. We find 
the following choice, with the MUB's arranged in a specific order, 
to be convenient ilL: 



A/o,Mi,--- ,M 



N, 



with 

M„ 
and 



{\m,r):r = 0,---,N-l} m = 0, • • • , TV - 1 

{|cx),r) : r = 0, • • • ,TV - 1} m = N 



|m, r) 




m = 0, • • • , AT - 1 

m — oo 



where V and W denote unitaries corresponding to gi and g2 re- 
spectively. For the sets obtained by replacing each vector in the 
ordered sets Mq, Mi, ■ ■ ■ , Mn by its projector we will use the no- 
tation 

Mti,Mi, ■ ■ ■ ,Mn- 
Explicit expressions for the projectors are given below [TT]: 

^^''^D{mj,j) m ^ oo. 

The MUB's at the extremes can clearly be identified as those asso- 
ciated with vertical and horizontal rays respectively. Again, as with 
striations, the vectors within each MUB, except for the last one, 
are related to each other (up to phases) by the action of Z)(1,0). 
Those in the last one are related to each other (up to phases) by 
D(0,1): 



\m, r){m, r\ 



£'(1,0): \m,r) 



D{0,1): \m,r) 




m ^ oo 
m — oo, 

\m, r — m), m ^ oo, 

I'm, r ~ 1) m — oo 



Quantum Nets: We have seen that each ray has a MUB associ- 
ated with it. We have also seen that the passage from one to 
another member of the striation by translation, has a complete 
correspondence with the passage from one vector in the MUB to 
another through appropriate displacement operators. We now as- 
sign a specific vector (more precisely the corresponding projector) 
in the MUB to each ray and assign other vectors in the MUB to 
the other parallel lines parallel to the ray using the correspon- 
dence mentioned above. This assignment of vectors to lines in the 
phase space consistent with translations leads to a Quantum Net. 
Clearly, as there are + 1 rays and for each ray N vectors to 



choose from the corresponding MUB's, there are iV^+^ different 
ways of assigning vectors to the hnes in phase space and hence 
]\fN+i (jigtinct quantum nets. Now given a quantum net we can 
generate another unitarily equivalent by applying one of the N'^ 
displacement operators to all the vectors. This then means that 
we can divide iV^+^ quantum nets into N'^~^ equivalence classes 
containing nets each. 

8. Phase Point Operators: Having assigned a vector to each line in 
the phase space we are now in a position to assign a Hermitian 
operator A{q,p) to each point {q,p) in the phase space by adding 
up all the projectors Px associated with the N + 1 lines A passing 
through that point and subtracting the Identity operator from the 
sum to make its trace equal 1: 

%P) = E^A-I; TrA{q,p) = l. 

A 

One can invert this relation to obtain the projectors associated 
with a line A in terms of the phase point operators 

Further using the mutually unbiasedness of the projectors one ob- 
tains: 

TT[A{q,p)A{q',p')]=NS,,g,Spy. 

From the above discussion it is clear that the totality {A} of iV^+^ 
phase point operators are obtained by picking one vector from each 
of the MUBs in all possible ways. 

9. Affine Planes: Given the phase point operator at say A(0,0) we 
can generate N'^ operators A{q,p) by applying all the N'^ displace- 
ment operators to it. Such a collection is taken to define an affine 
plane. We thus have the following correspondences between the 
structures at the phase space level and those at the level of the 
Hilbert space: 

Points {q,p) < — > Phase point operators A{q,p), 
Lines A •* — * Projectors, 
Phase Space < — > Affine Plane. 

The full set of N^~^ afHne planes can thus be obtained by dividing 
]\j-N+i jj^-j-Q subsets, related to each other by the Heisenberg- 
Wcyl group. 

10. Wigner Distributions: With the phase point operators at hand one 
can define Wigner distribution corresponding to a given density 
operator p as 

Wp{q,p)=Tr[pA{q,p)]. 
which, among others, have the desired marginals property: 

^ E W,{q,p)=Tv[pPx]. 
i.e. a Wigner distribution averaged over a line gives a probability. 
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3. A CONVENIENT LABELLING SCHEME FOR PHASE POINT OPERATORS AND 

AFFINE PLANES 

As discussed earlier, a selection of projectors, one from each A^o^-^ii ■ • ■ ,-Mn, 
defines a phase point operator A. We can therefore label a phase point operator 
uniquely by an array (rp, ri, • • • , rjv) : 

(ro,ri,-- - ,rjv). 

This permits us to think of each phase point operator as a vector in V/v+iC^at), 
an A^ + 1 dimensional vector space over Z^v- We now examine how they transform 
under the Clifford group. To this end we look at the action of the generators D{1, 0) 
and D(0, 1) of the Heisenberg-Weyl group as well as the action of the symplectic 
group. 

From (EB and ^ it follows that : 



i^(l,0) : (ro,ri,--- ,rA,)->(ro,ri,--- ,rjv) + (l,l,--- ,1,0), (29) 
i^(0,l) : (ro,ri,--- ,rA,)->(ro,ri,--- ,rjv)-(0,l,--- ,A-1,1). (30) 

Further, using and (PU]) it is straightforward (though somewhat tedious) to 
confirm that an arbitrary U{F) given by (|17p acts on |m,r)(m, r| as follows: 



U{F) : |m, r){m, r\ 



am-\-l3 



|oo, —77") (00, — 7r'| 



00, 5r) (00, 6r\ 



m ^ 00, 7m + 6^0 
m ^ 00, 7m + 6 = 
m = 00, 7 7^ 
m = 00, 7 = 0. 



On these translate into 

— (7TO — a)r f3~s-n 



m ^ N, 7m — a ^ 
m ^ N, 7TO — a = 
m = A, 7 7^ 
m = A, 7 = 



(31) 



(32) 



On the set of MUBs the action is that of a Mobius transformation. In fact the set 
of MUBs can be regarded as a projective line at infinity, added to the affine plane. 



These actions can be written out as matrices acting on r = (ro, 
for instance, for A = 3 we find for the generators of SL(2, Z3) 



where 



^{91) 







V 


r 


-^r'=Z^(.gi)r, 














w 


r 


r' = Uig2)r, 















1 






/o 










1 











; ^(92) = 








1 








1 











2 








Vo 








y 




V 











,r7v). Thus, 

(33) 
(34) 



(35) 



Note that these matrices have determinant equal to one and hence belong to the 
group SL(A + l,ZAr). 

Next we develop a similar notation for the affine planes. As discussed earlier, 
the set of phase point operators obtained by applying all elements of the Weyl- 
Heisenberg group to a fixed phase point operator defines an affine plane containing 
that phase point operator. This operation decomposes the set of A^+^ phase point 
operator into iV^^^ disjoint subsets — the affine planes, containing A^ points each. 
To develop a useful representation for the affine planes as cosets in Vn+iCZn) we 
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notice that the action of the Weyl-Heisenberg group on the phase point operators 
consists in adding vectors which are hnear combinations of cq and ei, where 

eo = (1,1,- •• ,1,0), (36) 
ei = (0,1,2,... ,7V-1,1). (37) 

Further, under the action of V and W these transform into each other: 

V : eo eo, ei ei - eo, (38) 
W : eo ^ ei, ei ~> -eo- (39) 

It therefore proves convenient to choose a basis which contains these two. Any 
choice would do. We choose the remaining ones as: 

et - (o,l^2^••• ,(7v-l)^o); fc = 2,... ,iv_i, 

eN = (0,0,0,--- ,0,1). (40) 

That this set of vectors is indeed a hnearly independent set can easily be checked 
using the properties of Van der Monde determinants. We shall denote the compo- 
nents of a vector (ro, ri, ■ ■ ■ , r^r) in the e-basis as [ao, ai, • • • , cxn] and the invertible 
matrix relating the two by S: 

r = 5a; 5= (ej,ef,--- ,e?;). (41) 

Under the action of the Weyl-Heisenberg group, the components in the e-basis have 
rather simple properties 

£'(1,0) : [ao,ai,--- ,aAr] ^ [ao-l-l,ai,- ,aAr], (42) 
D(0,1) : [Q;o,ai, • • • ,q;7v] [accki - 1, • • • ,aAr]. (43) 

Thus the collection [ao, cti, ■ ■ ■ , ajv], with ao and ai taking all values in Z^r with 
the rest fixed, defines an affine plane. We label it by the N — 2 coordinates 

["2, • • ■ , "AT-l, Oat]- 

The action of V and W is 

V: cx ^ a' ^ S-^U{gi)Sa, (44) 
W : cx^a' = S-^U{g2)SoL. (45) 

where the matrix 5* is given in (|4ip . 

If one is only interested in the actions of V and W on the affine planes then it 
suffices to look at the actions of the matrices S~^U{gi)S and S^^U{g2)S on the 
A'^ — 1 a's omitting ao and ai. These matrices have determinant equal to one and 
hence belong to the group SL(A — l,Zjv). 

To summarise, we have a representation of the SL(2,Z7v) action on the phase 
space through matrices belonging to SL(iV-|- 1, Z^r) ( SL(iV— 1, Z^v) ) acting on an 
A + 1 (A — 1) dimensional vector space over Zat whose elements are in one to one 
correspondence with the phase point operators (affine planes). A natural question 
to ask is whether or not one has an analogue of the symplectic form left invariant by 
SL(2,ZAr) action on the phase space. To pursue this question we look for bilinear 
quadratic forms of the type f7(r',r) = r'-^ Vl r, left invariant by the U{gi),U(g2) 
action on the r's, i.e we look for matrices 57 such that 

U^{gi)m{gi) = Q; U'^ {g2)m{g2) = n. (46) 

The antisymmetric matrix turns out to be 

n^(>;), (47) 
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where I denotes an TV-dimensional column with all entries equal to 1 and q is an 
N X N antisymmetric matrix with non-vanishing matrix elements 

q., = (48) 
In the e-basis, the symplectic form has a very simple structure: 

{0 I or j = or 1, 

J6^+JM 2<i,j<N-2, (49) 

-1 i = N -1 and j = N, 

which in turn suggest a basis in which Q has the canonical form. Thus, for instance, 
for iV = 5 the basis in which il has the canonical form 



IS 



I 

-I 



62,64, ^63,-65 ■ 



(50) 



Similarly for N — 7 we have 



1 1 

62,64,66, -65, -63,-67 . 
5 3 

4. Orbits under the Clifford Group 

Having learnt how the generators of the Clifford group act on phase point oper- 
ators and affine planes, our next task is to investigate how they arrange themselves 
into orbits under the action of this group. As is well known, to generate an orbit 
of a group G acting on a set we need to pick an element of the set and apply all 
elements of the group to it. We then pick another element in the set not contained 
in the earlier set and generate its orbit and so on until all elements of the set are 
exhausted. The quantities of interest are then (a) the number of orbits (b) the 
size of each orbit. The total number of orbits can be calculated using the Burnside 
Lemma: 

Number of orbits = — 0(g), (51) 

where \G\ denotes the order of the group and ^(g) the number of points of the set 
left fixed by the action of g. As g and the elements related to it by conjugation 
have the same fixed points, to calculate the number of orbits we need to know the 
number of elements in each conjugacy class and the number of points left fixed by 
a representative in each class. As to the size of the orbits, it is given by the ratio 
of the order of the group to the order of the stability group of the starting element 
of the set from which the orbit is built up by group action. Knowing the order of 
the subgroups, one can then deduce the possible sizes of the orbits. In summary, 
to answer the questions pertaining to the orbits, we need 

(1) knowledge of the conjugacy classes and the number of elements therein 

(2) expressions for the class representatives in terms of the generators 

(3) orders of the cyclic subgroups generated by of elements of a conjugacy class. 

In the following we collect together some relevant facts from the literature [T^ 
concerning these aspects for the group SL(2,Zjv) (also see the Appendix). 
In our discussion we will need to distinguish the sets 

Q = {x e : X = 2/^ for some y G Z^} (52) 

Q^{xeZ%:xiQ} (53) 
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where is the set of non-zero elements of Z^r (so Q is the set of quadratic residues 
and Q is the set of non-quadratic residues [H]). 

It is immediate that elements of SL(2,Z7v) having different trace must belong 
to different conjugacy classes. It turns out [12] that there is in fact exactly one 
conjugacy class for each value of the trace, except when the trace = ±2 in which 
case there are three. This gives us iV -f 4 conjugacy classes: the N "standard" 
classes 

with i== 0,1,... (TV 

C±2 - 



-1 

1 t 

1) and the four additional classes 
V ±2 



and 







[(o 





(54) 



(55) 



where v is any fixed element of Q. In these expressions \F\ denotes the conjugacy 
class containing F . For more details see Theorem [3] in the Appendix. 
For the standard classes we have 

{A^(A^-l-l) ift2-4GQ 
N{N-\) ift2_4GQ (56) 

\(N'^-\) if<2_4 = o 

where the notation \S\ means "number of elements in the set S"'. For the additional 
classes we have |C'±2| = (^^ — l)/2 and |D±2| = 1- It may be worth noting that 
the classes Ct for which <^ — 4 e Q are precisely the ones whose elements are 
diagonalizable (apart from the two classes D±2)- 

For trace = 0, ±1,±2 it is easy to calculate the orders of the cyclic subgroups 
generated by the elements in each conjugacy class for arbitrary iV. They are listed in 
Table [T] (to derive the result for elements of C±2, C±2 consider matrices of the form 



class 


Co 


Ci 




C2 


C2 


D2 




C-2 


D-2 


order 


4 


6 


3 


N 


N 


1 


2N 


2N 


2 



Table 1 . The order of the cyclic subgroups generated by elements 
of some special conjugacy classes, for all N. 



{^0 ±i))- ^^'^ remaining conjugacy classes the orders of the cyclic subgroups 
depend on number theoretical details. We have worked out the orders for N < 19. 
They are listed in Table [21 Since we know how many elements there are in each 





C3 


C4 


C5 


Ce 


C7 


Cs 


C9 


C'lo 




C12 


Cl3 


Ci4 


Ci5 


C16 


TV = 7 


8 


8 


























N =11 


5 


10 


12 


12 


5 


10 


















TV = 13 


14 


12 


14 


14 


7 


7 


12 


7 














TV = 17 


18 


18 


16 


8 


9 


16 


16 


18 


8 


16 


9 


9 






iV= 19 


9 


5 


10 


20 


9 


20 


9 


18 


20 


18 


20 


5 


10 


18 



Table 2. The order of the cyclic subgroups generated by elements 
of all the conjugacy classes, for TV < 19. 



conjugacy class, the number of distinct cyclic subgroups can easily be worked out. 
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order of 












number 


cyclic 


fixed 










of 


sub- 


points in 




fixed 


N 


Class 


elements 


group 


plane 


fixed As 


planes 


3 


Co 


6 


4 


1 


1 


1 




Ci 


4 


6 


1 


1 


1 




Ci 


4 


6 


1 


1 


1 




Di 


1 


2 


1 


1 


1 




C2 


4 


3 


3 


9 


3 




C2 


4 


3 


3 


9 


3 




D2 


1 


1 


9 


81 


9 


5 


Co 


30 


4 


1 


1 


1 




Ci 


20 


6 


1 


1 


1 




C2 


12 


5 


5 


25 


5 




C2 


12 


5 


5 


25 


5 




D2 


1 


1 


25 


15625 


625 




C3 


12 


10 


1 


1 


1 




C3 


12 


10 


1 


1 


1 




D3 


1 


2 


1 


1 


1 




Ci 


20 


3 


1 


25 


25 


7 


Co 


42 


4 


1 


1 


1 




Ci 


56 


6 


1 


1 


1 




C2 


24 


7 


7 


72 


7 




C2 


24 


7 


7 


72 


7 




D2 


1 


1 


72 


78 


76 




C3 


42 


8 


1 


1 


1 




Ci 


42 


8 


1 


1 


1 




Cb 


24 


14 


1 


1 


1 




c^ 


24 


14 


1 


1 


1 






1 


2 


1 


1 


1 




Ce 


56 


3 


1 


72 


72 



Table 3. Fixed points under the action of the elements in 
SL{2, Zjy), N < 7. The action can be on an affine plane, on 
the set of phase point operators, and on the set of affine planes. 



In particular we find that there is a unique subgroup of order 2, generated by the 
element — / £ D^2\ N + 1 (Sylow) subgroups of order TV, generated by elements of 
C2, C2; and + 1 subgroups of order 2N generated by elements of C_2, C'-2- 

We can now compute quantities of interest concerning the orbits generated by 
Clifford group action on phase point operators and afhne planes. Using the results 
in the previous section we have done this for N = 3,5,7. The results are listed 
in Table [H Some things can be said in general. Notably cyclic subgroups of even 
order have only one fixed point. It appears that cyclic subgroups of an odd order 
m always leave exactly affine planes fixed, where + 1 = km + m' , m' < m. 
(We proved this for subgroups of order 3, and checked it in all cases for N < 19.) 

On the afhne plane the number of orbits is always 2. On the set of affine planes 
the number of orbits can be found using Burnside's lemma. The results for iV < 11 
are listed in Table H) They should be compared with the estimate N^~^/N^ 
which for N — 5,7, 11 gives 5,343, 19487171 respectively. The estimate is rather 
good. For = 3 the sizes of the orbits are 1 + 8 = 32, and for N = 5 they are 
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N 


3 


5 


7 


11 


No. of orbits 


2 


11 


360 


19650810 



Table 4. Number of orbits of SL{2, Zm) when acting on the set 
of afhne planes. 



1 + 24 + 40 + 40 + 40 + 40 + 40 + 40 + 120 + 120 + 120 = b"^. For iV = 7 the results 
are too numerous to be given here. As expected, each term in the sum is a divisor 
of the order N{N'^ — 1) of the group. We observe that for all N there is a unique 
singlet, and a unique {N"^ — l)-plet consisting of afRne planes left invariant by some 
Sylow subgroup.. 

5. Spectra of Phase Point Operators, Complex Conjugation and the 

Extended Clifford Group 

We now turn to the spectra of the phase point operators. Obviously those in the 
same affine plane have the same spectra and hence we need to focus only on the 
affine planes. For N = 3 one finds two distinct spectra: (1, 1, —1) and ($, 1 — $, 0) 
where $ = (1 + \/5)/2 = 1.61803.., the Golden Ratio. No surprises here. 

For N = 5 one finds only 9 distinct spectra as against the naive expectation of 
11, as listed in Table [S] For N = 7 one gets 210 distinct spectra (too numerous 



spectrum 


no. of occurrences 


{-l.,-l.,l.,l.,l.} 


1 


{-1., -0.61803, 0,1., 1.61803} 


24 


{-0.94658, -0.5169, -0.18438, 0.93842, 1.70944} 


120 


{-0.90932, -0.48701, 0, 0.46853, 1.9278} 


120 


{-0.90039, -0.64018, -0.14531, 1.06785, 1.61803} 


40 


{-0.83726, -0.58152, -0.09576, 0.6287, 1.88584} 


120 


{-0.83607, -0.81, 0, 1.05469, 1.59139} 


80 


{-0.79859, -0.36221, 0, 0.10661, 2.05419} 


80 


{-0.70281, -0.61803, -0.13294, 0.48666, 1.96712} 


40 



Table 5 . Possible spectra of the phase point operators in an affine 
plane, TV = 5. 



to reproduce here). In both cases one finds that sizes of the spectral orbits do not 
divide the order of the group but rather twice the order. This leads one to suspect 
that as far as the spectra are concerned, the relevant group is not the Clifford group 
but rather an augmented Clifford group. The extra symmetry is provided by the 
operation of complex conjugation which also maps the MUB's into themselves. The 
Extended Clifford Group [5] includes the complex conjugation operation and has 
twice as many elements as the Clifford group. This is consistent with the sizes 
of the spectral orbits found. With this in mind we investigate and encode the 
action of complex conjugation on the MUB's, so as to have explicit actions of all 
the operations on the MUB's relevant for discussing the spectra of the phase point 
operators which in turn enable us to classify distinct Wigner distributions. 



12 



Complex conjugation, denoted by C, is represented by an anti-unitary operator. 
Under C the MUB's transform into each other as follows: 



which implies 



C : \m, r) 



C: 




m 7^ oo 
m — oo 



m = 0, 1,...,7V 

m = TV 



and 



0,1, 



(57) 



(58) 



(59) 



These actions, when translated into actions on the a's, have the pleasant feature 
that under complex conjugation the a's alternately change signs: 



C : [ao, ai, • • • , un] [ao, -ai, ^2, 



-oln\ 



(60) 



Thus complex conjugation acts on the phase point operators [q;o,q;i, • • • , aAr] by 
a diagonal matrix diag(l, —1, 1 • • • , —1), with determinant — 1 if = 4fc + 1 and 
1 if iV = 4fc + 3 . On the afhne planes [012, as, • • ■ tChn], the action is again by a 
diagonal matrix but with determinant 1 if iV = 4fc + 1 and —1 if = 4fc + 3. 

As described in ref. [S] the extended Clifford group is obtained by admitting 
matrices in ESL(2,ZAr) (the group consisting of 2 x 2 matrices with entries in Z^r 
and determinant = ±1). The conjugacy classes of ESL(2, Z^r) are described in the 
Appendix to this paper. Unlike SL(2,ZAr) there is a distinction between the cases 
A^ = 1 (mod 4) and A^ = 3 (mod 4) . It turns out that for all values of A^ we have 
the 2N "standard" classes 

' ^0 -A' 
1 t 



Ca.. 



(61) 



where A = ±1 and t = 0, 1, 



. , (A^ — 1), and the 2 additional classes 

D±2 = 







[(0 


°01 



(62) 



If A^ = 3 (mod 4) these are the only conjugacy classes (see Theorem [5] in the 
Appendix). If, however, A^ = 1 (mod 4) there are another 6 conjugacy classes, 
which we denote C±2, C±2i, D±2i (see Theorem [T] in the Appendix). 

All elements of ESL with determinant —1 are represented on the symplectic 
vector space of affine planes, Vn-iC^n), by matrices A that obey 



A^flA = -n . (63) 

They are therefore anti-canonical transformations. In Table [H] we list the number 
of fixed points for the conjugacy classes of ESL(2,ZAr) with determinant —1. For 
all A^ one can show that elements of order 2 leave A^(^~i)/2 affine planes fixed. 
Because the transformations are anti-canonical it follows that these fixed planes 
form Lagrangian subspaces of Vn^iC^n)- Furthermore we believe that if an element 
with determinant —1 squares to an element leaving A^^*"' As fixed, then that element 
itself leaves A^*^ As fixed; we checked this for A^ < 11. 

Again we can use Burnside's lemma to compute the number of orbits under 
ESL(2, Zat). The results for A^ < 11 are hsted in Table H For large A^ the number 
approaches the estimate N'^~-^ /2N^ . For N < 7 we know the number of distinct 
spectra, and we find complete agreement with the number of orbits. 
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order of 












number 


cyclic 


fixed 










of 


sub- 


points in 




fixed 


N 


Class 


elements 


group 


plane 


fixed 


planes 


3 


C'-i.o 


12 


2 


3 


3^ 


3 






6 


8 


1 


1 


1 






6 


8 


1 


1 


1 


5 


C-1.0 


30 


2 




5^ 


5^ 




C-1.1 


12 


20 


1 


1 


1 




Ci 


12 


20 




1 


1 




Di 


1 


4 




1 


1 






20 


20 


1 


1 


1 






20 


12 


1 


1 


1 




C-iA 


12 


20 


1 


1 


1 






12 


20 


1 


1 


1 




Di 


1 


4 




1 


1 


7 




56 


2 




3^1 


—n 

73 




C-is 


42 


16 




1 


1 




C-1^2 


56 


6 




7 


7 






42 


16 




1 


1 




C-iA 


42 


16 




1 


1 






56 


6 




7 


7 






42 


16 




1 


1 



Table 6. Fixed points under the action of elements in ESL{2, Zn) 
with determinant —1,N<7. 



N 


3 


5 


7 


11 


No. of orbits 


2 


9 


210 


9833460 



Table 7. Number of orbits of ESL{2, Z^) when acting on the set 
of aflane planes. 



6. Summary 

To summarise: 

• For A'' an odd prime, a compact, mathematica friendly way of representing 
phase point operators and affine planes as vectors and cosets in an A'' + 1 
dimensional vector space over Z^r is given. 

• In this representation the SL(2, Z^) actions on phase point operators ac- 
tions by matrices in SL(A'' + 1, Zjv). 

• A convenient basis, factoring out the Heisenberg-Weyl group actions, is 
introduced in which SL(2, Z^) actions on the affine planes can be described 
in terms of matrices in SL(A'^ — 1, Zj^). 

• The action of complex conjugation on the phase point operators is by ma- 
trices in GL(A'^ -|- l,Zisr) with determinant —1 if N = 4fc -|- 1 and 1 if 
A^ = 4fc + 3. 

• The action of complex conjugation on the affine planes is by matrices in 
GL(A'' - 1, Zn) with determinant= 1 if A/' = 4fc -|- 1 and -1 if A'' = 4fc 3. 
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• The {N — l)-dimensional vector space of afSne planes is symplectic. El- 
ements of ESL(2, Zn) with determinant 1 preserve the symplectic form, 
elements with determinant —1 are anti-canonical. 

• For N < 7 the number of distinct spectra in the phase point operators 
equals the number of orbits under the group ESL(2,ZAr). 

The work reported here suggests a number of avenues for future research. In the 
first place one would like to know if the degeneracies of the phase point operators 
are completely determined by ESL(2,ZAr) for every odd prime dimension N. Be- 
yond that one would like to know if the result generalizes to the case of arbitrary 
odd prime power dimension. Finally, there is the question of even prime power di- 
mensions. In that connection let us observe that Gibbons et al [5j find 4 similarity 
classes when iV = 4. However, when we computed the spectra we found only the 
3 distinct spectra tabulated in Table [HI The cause of this additional degeneracy 



spectrum 


no. of occurrences 


{-0.896802, -0.14204, 0.278768, 1.76007} 


384 


{-0.866025, -0.5, 0.866025, 1.5} 


320 


{-0.5, -0.5, 0.133975, 1.86602} 


320 



Table 8. Possible spectra of the phase point operators in an affine 
plane, iV = 4. 



cannot be the same as the cause in odd prime dimension since when N is an even 
prime power the fact that —1 = 1 (mod 2) means that SL(2,FAr) = ESL(2,F7v). 
There must therefore be some other explanation, which it would be interesting to 
investigate. 
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Appendix: Conjugacy Classes of ESL(2,FAr) 

In this appendix we deduce the conjugacy classes for the group ESL(2, Fjv) where 
A — p' is a power of an odd prime number p and is the Galois field [14j having 
A elements (in the main text we only need the result for the case ^ = 1, F^r = Z^r; 
however, the result for arbitrary I is no more difficult, and it may be useful for 
future developments). For the sake of completeness, and the convenience of the 
reader, we also describe the conjugacy classes of the group SL(2,Fjv) (also given in 
ref. dH 

Let F^ be the set of non-zero elements of Fat, and define 

Q ^ {x e : X ^ y"^ for some y G ¥*j^} (64) 
Q^{xe¥%:xiQ} (65) 
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Let 6 he a primitive element [Tl] for ESL(2, F^r). Then it is easily seen that 9" G Q 
if and only if s is even. Consequently Q and Q each contain exactly {N — l)/2 
elements. Also the fact that ^(^^i)/^ = — i means that — 1 G Q if and only if 
N = 1 (mod 4). 

We begin by proving three preliminary results. 



Lemma 1. For all x e 



jV~l 
4 



if N = I ( mod 4 ) and x £ Q 



\Qr\{Q-x)\^ ifN=l(modA)andxGQ (66) 

^ if N = 'd (mod 4 ; 

where the notation \S\ means "number of elements in the set S". 

Proof. Let x be any fixed element E F^. Define 

5, = (xQ) n{Q~-x) (67) 

= (xQ) n{Q-x) (68) 
and let fx ■ F|^ F^ be the map defined by 

fxiy)^x^y~\ (69) 

It is easily verified that fx is a bijection and (using the fact that QQ = QQ = Q 
and QQ = Q) 

f.{Sx)^Tx fx{Tx) = Sx. (70) 

Since Sx, Tx are bijective images of each other we deduce that \Sx\ — \Tx\- 

Suppose, now, that N = \ (mod 4). Then —1 € Q which means —x ^ xQ for all 
non-zero x. Consequently \Sx U Tx\ — \xQ\ = {N — I)/2, implying 

1^.1 = 1^.1 = ^ (71) 

Suppose, on the other hand, that = 3 (mod 4). Then — 1 G Q which means 
—x G xQ for all non-zero x. Consequently jS'^rUTj;! — \xQ\ — 1 = {N ~ 3)/2, 
implying 

- 3 

l^.l = \Tx\ = (72) 

To complete the proof suppose, first of all, that x E Q. Then xQ = Q, implying 

I - / - M , , f ^ N = 1 (mod 4) , , 

\Qn(Q-x)\= Sx^ I ^ 73 
^' ' ' A = 3(mod4) ^ ^ 

Suppose, on the other hand, that x E Q. Then xQ = Q, implying 

A — 3 

\Qn {Q-x)\ = — \Qn{Q-x)\ 



A^ = 1 (mod 4) 
A^ = 3 (mod 4) 



(74) 
□ 



Lemma 2. For all fi G F^ and v G F^ there exists q E Q U {0} such that 

W + G Q U {0} (75) 
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Proof. Case 1: ^ € Q. We have |(5U{0} + i^| > \Q\. Consequently the set 
{Q U {0} + v) r\{Q\J {0}) is non-empty. We can therefore choose q' & Q \J {0} 
such that + I' e Q U {0}. Then q = fi^^q' has the stated property. 

Case 2: ^ ^ Q. If = we can choose q = Q. If, on the other hand, v ^ we 
note that it follows from Lemma [T] that | (0 + ^) 0| < lOl- So there exists q' £ Q 
such that q' + E QU {0}. Then q = n^^q' has the stated property. □ 



F=\Z (76) 



Lemma 3. Let 

w- i 

he any matrix G ESL(2,FAr). Let A = deti^ and t = tr F. Then 

(1) - 4A 7^ there exists S G ESL(2,FAr) such that 

F^s(^^ -f^y-' (77) 

(2) Ift^~AA = andP^O then, for any u e there exists S G ESL(2, Fjv) 
such that 

^ = t)*" 

(3) //t^ - 4A = BTirf 7 7^ then, for any u G F^, there exists S G ESL(2, Fa,) 
such that 

F^S^^l (79) 
Moreover in every case the matrix S can be chosen so that det S — A. 
Proof. Let k G F^ be arbitrary and look for a matrix 

S=(^ (80) 

such that 

F^S(1 -t)^- (81) 
It is easily seen that S has this property if and only if 



which means S must be of the form 



^ ax+l3y 



^-[y l^Sy) (83) 

Taking into account the requirement det — A the problem thus reduces to the 
problem of finding x, y such that 

7^^ + {5 - a)xy - = fcA (84) 
Case 1: t^ — 4 A ^ 0. Setting = 1 in Eq. ((5i| the condition becomes 

7a;2 + {5- a)xy ^ f3y^ = A (85) 

If 7 = this equation has the solution x — {A + f3) / {d — a), y = 1 (note that the 
fact that t'^ — 4A means 6~a ^ 0). If, on the other hand, 7 7^ we can rewrite 
the equation in the form 

^^(^y^(tl^4A)^^±47A 
27 / 472 

It follows from Lcmma[5]that there exist r, s G Fjv such that {t'^ — 4:A)s'^+4^A = r^. 
The equation then has the solution x — (r ^ s{d — a)) /(2"f), y = s. 
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Case 2: - 4A = and (3 0. Setting k = -(iu^ in Eq. (HH) the condition 
becomes, after rearranging, 

2 ,2+2 



y 



{S — a)x 
2/3 



which has the sohition x = 1, y = (ut(3 + {5 — a)) / (2/3). 

Case 3: — 4 A ~ and 7 7^ 0. Proved in the same way as Case 2. 



(87) 



□ 



We are now ready to deduce the conjugacy classes. We write F ~ G if F is 
conjugate to G, and use the symbol [F] to denote the conjugacy class containing 
F. Define 

" ^0 -A' 
1 t 



where A = ±1 and t is any element of . 
conjugacy classes. We also need to consider the classes 



''n- We refer to these as the "standard" 











and 




t J. 







Ml 




t/2j 



(89) 



where is any fixed element of Q and t is such that = ±4. We then have 



Theorem 1. Let N — 1 (mod 4). Let i be one of the two elements o/F^v with the 
property — —1. Then the conjugacy classes o/ESL(2,Fjv) comprise 

(1) The 2N standard classes CA,t with A = ±f and t G F^r. 

(2) The 8 additional classes C±2, C±2i, D±2, L)±2i- 



Let F 



(;f) 



(1) Ift^-iA^O 

(2) //t2-4A = 



and let A = detF, t = tiF. Then 
F e GA,t 



(90) 



Fe 



(note that if t^ - 4 A 
/? G Q and 7 G Qj. 



Ga,* if(3or-feQ 

Ct ifPor-feQ (91) 

Dt z//3 = 7 = 

it cannot happen that /? G Q and 7 G or that 



Remark. So the determinant and trace completely fix the conjugacy class except 
when - 4A = 0. 

Proof. To prove that the classes listed are disjoint note that F ^ G implies det F = 
detG and trF = trG, and that the classes Dt each consist of a single element. It 
is therefore enough to show that 



GA.t nCt = 



(92) 



when - 4A = 0. To see this suppose F = (7 a) ^ C'i,±2- Then there exists 

(93) 



S={l^)e ESL(2, Fjv) such that 

F^S 



-1 

1 ±2 



Performing the algebra we find 



a P 
7 S 



As(wT2/)^ 



-AsizTx){wTy) ± 1, 



(94) 
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where As = det S. It follows that 7 G QU{0}. We conclude that ° ^22") C'i,±2- 

The proof that the classes C_i^±2i and C±2i are disjoint is similar. 

The remainder of the theorem is a straightforward consequence of Lemma [31 

□ 

If = 3 (mod 4) there are fewer conjugacy classes. The reason is that in 
this case — 1 G Q. Consequently the classes D±2i and C±2i do not exist. Also 
E Q, implying —v — for some q € F^. So the similarity transformation with 

S = (0 -i/g) takes ( " "±^" ) to ( ? ±1 ) , implying C±2 = Ci^±2- We thus have 

Theorem 2. Let iV = 3 (mod A). Then the conjugacy classes o/ESL(2,Fjv) com- 
prise 

(1) The 2N standard classes C'^^t with A = ±1 and t G F^r. 

(2) The 2 additional classes D±2. 

Let F be any element o/ESL(2,Fjv) and let A — det_F, t = tr_F. Then 

(1) IfF^±I 

F e Ca,* (95) 

(2) IfF = ±I 

F e £'±2 (96) 

Remark. Here too the determinant and trace completely fix the conjugacy class 
except when t'^ — 4 A = 0. 

Finally, for the sake of completeness, let us note that Lemma [3] can also be used 
to deduce the conjugacy classes of SL(2,F7v) (which are given by, for example, 
Gehles [H]). Since A = 1 for all elements of SL(2,Fjv) we use the simplified 
notation (for the standard classes) 

1 t 



Ct 



(97) 



and (for the additional classes) 

L^±2 = 



and 



D±2 = 







[(0 


°01 



(98) 



where v is any fixed element of Q. We then have 
Theorem 3. The conjugacy classes o/SL(2,FAr) comprise 

(1) The N standard classes Ct, for arbitrary t G Fjy. 

(2) The 4 additional classes D±2, C±2- 

LetF^(^"1^€ SL(2, Fat), and let t = iiF. Then 
(1) Ift^±2 



FeCt 



(99) 



(2) Ift = ±2 



Fe 



if —(3 or J G Q 

if -(3 or J eQ (100) 
13 G Q and ^ G Q, or that 



(note that if t — ±2 it cannot happen that 
-(3 £Q and 7 G Qj. 

Remark. So the trace completely determines the conjugacy class except when t = 
±2. Note that for the group SL(2,Fjv) there is no distinction between the cases 
iV = 1 (mod 4) and iV = 3 (mod 4). 

Proof. Similar to the proof of Theorem [T] □ 
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